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Abstract 

We investigate the physical mechanism for color production by isotropic nanostructures with 
short-range order in bird feather barbs. While the primary peak in optical scattering spectra results 
from constructive interference of singly-scattered light, many species exhibit secondary peaks with 
distinct characteristic. Our experimental and numerical studies show that these secondary peaks 
result from double scattering of light by the correlated structures. Without an analog in periodic 
or random structures, such a phenomenon is unique for short-range ordered structures, and has 
been widely used by nature for non-iridescent structural coloration. 
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Structural colors have attracted much attention in a wide variety of disciplines {l| . They 
originate from physical interactions of light with nanostructures. Recent studies have focused 
on ordered structures in the natural world which produce iridescent colors 2, 3,0]. However, 
nature has also used extensively quasi-ordered structures to create vivid colors that are 
weakly iridescent, e.g., the spongy nanostructures in the feather barbs of numerous birds 
[5I. In the 19th century, lack of pronounced iridescence led to the hypothesis that the colors 
were produced by wavelength selective scattering from single particles. Raman was the 
first to cast doubt on this hypothesis |6|, which was later falsified by Dyck Prum et 



al. proved the existence of structura 
medullary keratin in feather barbs 



order by Fourier analysis of electron micrographs of 
]. In the last ten years, the previously unappreciated 
class of quasi-ordered nanostructures have been identified also in bird skins, mammal skins, 
dragonfly cuticles, and butterfly scales HQ. 

Although it is now recognized that light scattering is organized to contribute effectively to 
coloration of quasi-ordered structures, there are significant spectral features whose origin re- 



mains elusive 



quasi 

m 
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In previous studies, optical reflectance spectra are predicted from 



the structural Fourier spectra with the assumption of single scattering 1^. The dominant 
length scale of structural correlation gives one peak in the predicted reflectance spectrum. 
However, many species we have measured feature two reflection peaks, in agreement with 
prior observations . [1^. What is the origin of the second peak? Does multiple scattering 
play a significant role in coloration? How does the coloration of quasi-ordered structures 
differ from that of ordered structures? Answering these questions is essential to identify 
the physical mechanism for coloration of the quasi-ordered structures that are widely dis- 
tributed in nature and clearly function in animal communication. Moreover, understanding 
the effects of short-range structural order on multiple scattering of light is fundamental to 
the physics of light interaction with complex media. In practical terms, the wavelength 
dependence of light scattering from quasi-ordered structures has not been taken advantage 
of technologically, and may have useful advantages for photonic coatings across many ar- 
eas such as the production of color for paint, cosmetics and textile industries and photon 
management in solar cells. 

In this Letter, we investigate the physical mechanism for coloration of nanostructures 
with short-range order in bird feather barbs by performing angle- and polarization-resolved 
scattering spectrometry. The primary peak in optical scattering spectra maintains the po- 
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larization of incident light and its spectral shift with observation angle agrees with single 
scattering predictions. The short-range structural order leads to phase correlation of scat- 
tered waves and their constructive interference produces a narrow spectral peak. For many 
avian species, the scattering spectra have a secondary peak which is depolarized and has 
different angular dispersion from the primary peak. We show that the secondary peak is 
caused by interference of light scattered twice by the nanostructures. This result is unex- 
pected from the general belief that multiple scattering of light simply broadens the single 
scattering peak. There is no analogous phenomenon in periodic structures where multiple 
diffraction or high-order diffraction does not create additional peaks. 

The color-producing quasi-ordered nanostructures within the medullary cells of bird 
feather barb rami belong to two morphological classes [s], [isl. Channel-type nanostruc- 
tures consist of /5-keratin bars and air channels in elongate and tortuous forms. Sphere-type 
nanostructures consist of spherical air cavities in a /?-keratin matrix. Our recent studies 



suggest these nanostructures are se^ 
tein from the cytoplasm of the cell 



f-assembled during phase separation of /?-keratin pro- 



In this work we have performed small angle X-ray 



scattering (SAXS) and optical measurements on both types of structures from many bird 
species. Since the results are similar, we present here the data from turquoise wing feathers 
of Comcias benghalensis (C. benghalensis) , the same bird species that Raman studied over 
seventy years ago [Fig. [U^a)]. 

(a) (b) 




FIG. 1: (a) Coracias benghalensis (Indian Roller) whose wing feather color was studied by Raman 
in the 1930s. Photo credit to Pascal Rapin (reproduced with permission) (b) Scanning electron 
micrograph of the feather barbs showing channel-type nanostructures with /J-keratin rods and air 
channels. The scale bar = 500nm. Inset: small-angle X-ray scattering (SAXS) pattern of the 
feather barbs that exhibits a ring. The scale bar = 0.05 nm~^ 



The barb nanostructures of C. benghalensis are a 3D network of keratin rods of ~100 
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nm diameter and air channels of similar width [Fig. [U^b)]. The SAXS pattern exhibits a 
ring [inset of Fig. [Hb)], implying the nanostructure is isotropic. Exploiting the rotational 
symmetry, we average the scattering intensity azimuthally to obtain the intensity J as a 
function of the magnitude of spatial vector q [Fig. [2|^a)]. The strong peak of I{q) at gg 
= 0.028 nm^^ reveals the existence of a dominant length scale s for structural correlation. 
s = 27r/go = 225 nm is approximately half of the wavelength of blue light. The full width 
at half maximum (FWHM) Aq of the peak in I{q) reflects the range C, of spatial order, 
^ = 271/ Aq = 1.16 fim. Since ^ is only 5 times of the spatial period s, the order is short- 
ranged. Such structures can be called amorphous photonic materials, in analog to amorphous 
electronic materials with only short-range order. 

(a) (b) 
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FIG. 2: (a) Azimuthal average of SAXS intensity (blue dashed line) in comparison to nearly 
backscattered light intensity at normal incidence (black solid line) plotted as a function of the 
magnitude of spatial vector q. The optical scattering spectrum has a secondary peak, (b) Scattered 
light intensity of co-polarization (black solid line) and cross-polarization (blue dashed line) as a 
function of wavelength. The linearly polarized white light is incident at normal, (p = 0° . The 
detector collects the scattered light at an angle of 9 = 10° from the surface normal. The polarization 
direction of the incident light is parallel to the scattering plane that is formed by the surface normal 
and detector. Similar results are obtained when the polarization of incident light is perpendicular 
to the scattering plane. 

We have measured the scattered light from the feather barbs as function of wavelength, 
polarization, sample orientation, incident angle and observation angle. The cylindrical- 
shaped feather barbs are mounted at the rotation center of a goniometer, with their axes 
perpendicular to the rotation axis of the goniometer. Collimated white light is incident on 
the sample at an angle from the surface normal. Scattered light is collected by a lens and 
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focused to an optical fiber that is connected to a spectrometer. The measured spectra of 
scattered hght are normahzed by the spectrum of incident hght. The angular resolution of 
detection is about 5°. Linear polarizers select the scattered light whose polarization is either 
parallel (co-polarization) or perpendicular (cross-polarization) to that of incident light. 

To further confirm the origins of the two peaks, we have measured their angular depen- 
dences in a series of experiments. In the first experiment, the direction of incident light 
and the detector position are fixed. As the sample rotates, the frequencies of the primary 
peak and secondary peak in the scattering spectra do not change. Their independence of 
the sample orientation is consistent with the isotropic nature of the nanostructures. Next 
we fix the sample orientation and the incident angle of white light, and move the detector to 
measure the spectra of light scattered into different directions. Both the primary peak and 
secondary peak shift in wavelength with the detection angle, but in the opposite directions. 
We also change the incident angle and repeat the measurement. Analysis of all the data 
reveals that the frequencies of primary peak and secondary peak do not depend separately 
on the incident angle or the observation angle, but on their difference, that is, the angle 
9 between the directions of illumination and observation [inset of Fig. [3]^a)]. As shown in 
Fig. [3t^a), the primary peak shifts to shorter wavelength and the secondary peak to longer 
wavelength when 9 increases. We have also measured the reflection spectra of white light for 
different incident angle (j). A comparison to the scattering spectra reveals that the reflection 
peaks at an incident angle have the same frequencies as the scattering peaks at angle 
9 = 2(f) [Fig. [3](a)]. Since the angle between the incident beam and the reflected beam is 
9 = 2(f), their spectral coincidence implies the scattering peaks and reflection peaks have the 
same physical origin. Therefore, specular reflection can be regarded as a special case of light 
scattering from the quasi-ordered structures. 

For single scattering of light from the incident wavevector kj to a wavevector ko, ko — kj = 
q, where q is a spatial vector of the structure. 9m is the angle between — kj and ko [inset of 
Fig. [3] (b)]. The scattering is elastic, |kj| = |ko| = k. According to the inset diagram of Fig. 
[3] (b), the magnitude of q is 



where k = 27me/X, is the effective index of refraction, and A is the vacuum wavelength. 
We compute the q values for the primary peak and secondary peak in the measured scat- 
tering spectra. Light refraction at sample surface is taken into account by converting the 
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FIG. 3: (a) Inset: a schematic showing the white Ught incident on the feather barb at an angle (p 
from the surface normal is scattered or reflected to a direction of angle 6 from the incident beam. 
In the specular reflection measurement, 9 = 2(j). Main panel: Wavelengths of the primary (squares) 
and secondary (circles) peaks in optical scattering spectra (solid symbols) and specular reflection 
spectra (open symbols) as a function of 9. The top horizontal axis represents the incident angle 
4> in the specular reflection measurement, (b) q values for the primary peak (squares, solid line) 
and secondary peak (circles, dashed line) as a function of scattering angle 9. The discrete symbols 
represent experimental data and the lines calculations. Inset: single scattering diagram. 

experimental angle 9 outside the sample to 9m inside the sample [l^. As shown in Fig. [3] 
(b), the q value for the primary peak of light scattering is the same for different 9. It is 
equal to go of the SAXS data [Fig. EJ^a)] when rie = 1.25. This value corresponds to a vol- 
ume fraction of air of 57% in the keratin matrix, which matches the estimate from electron 
micrographs. Therefore, the primary peak in the optical scattering spectra is produced by 
single scattering process involving the dominant spatial vectors of the nanostructure. 

Unlike the primary peak, the q value for the secondary peak changes with 9. We also 
compare the optical scattering spectra to X-ray scattering results by plotting the nearly 
backscattered light intensity {9 = 10°) against q in Fig. EJ^a). While the primary peak of 
light scattering matches that of SAXS, the secondary peak finds no correspondence in the 
SAXS data. This discrepancy is attributed to stronger scattering of light than X-ray by 
the nanostructures, which results in multiple scattering of light. The lowest order multiple 
scattering is double scattering, i.e., the incident light is scattered successively twice before 
exiting the sample. As shown schematically in Fig. H] (a), light is scattered from kj to k^, 
and then from k^ to k^. The two successive single scattering events are considered to be 
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FIG. 4: (a) Two representative diagrams of double scattering. Light is scattered successively from 
kj to km, and then to kg. qi = k^ — kj and q2 = kg — k^ are spatial vectors of the structure. 
The angle Om between — kj and ko is the same in the two diagrams. The single scattering angles 
ai and a2 are smaller than 90° in the upper diagram and larger than 90° in the lower diagram, 
(b) Calculated double scattering spectra for different scattering angles 9 when the incident angle 
(f) is fixed at 0°. The values of 9 are marked next to the curves. 



uncorrelated. k.^ — kj = qi, and ko — k^ = q2, where qi and q2 are spatial vectors of the 
structure. Light scattering is elastic, |kj| = \km\ = |ko| = k. When the directions of kj 
and ko are fixed, k^ can be in multiple directions, as illustrated in Figs. ID^a). ai and a2 
denote the angles for two successive single scattering processes. As ai and a2 increase, k 
must decrease in order to keep qi and q2 close to go- 

We calculate the double scattering spectra of light from the Fourier power spectrum 
/(q) of the structure, /(q), which is measured directly by SAXS, tells the spatial vectors 
existing in the structure and their strength. Since the structure is isotropic, /(q) = I{q). 
For simplicity, we ignore the polarization of light. The probability of incident light with kj 
being scattered to k^ is proportional to /(qi = km — kj), and the probability of successive 
scattering from km to ko is proportional to /(q2 = ko — km). The summation of /(qi)/(q2) 
over all possible directions of km gives the probability of incident light being scattered 
to direction 6m- In the calculation we sample 1000 directions of km which are uniformly 



distributed over in solid angle 



15j. For comparison with the experimental data, we take 



into account light refraction at sample surface and sum the scattering spectra over a 5° 
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collection angle. 

Figure Hl^b) shows the calculated double scattering spectra when the incident light is 
normal to the sample surface. As 9 increases, the double scattering peak is broadened and 
split into two subpeaks. The subpeak at longer wavelength is stronger, and it red-shifts with 
increasing 9. The broadening of double scattering spectrum can be understood qualitatively 
by approximating I{q) as a delta function 5{q — go)- When 9m = 0, ai and 02 can only be 
90°. Only one wavelength is allowed and its value can be obtained by replacing 9m in Eq. [1] 
by «! or a2- Once 9m > 0, ai and 02 can take multiple values. The wavelength is no longer 
unique. The larger the 9m, the wider the range of possible values for ai and 02, and the 
broader the double scattering spectrum. 

Experimentally the efficiency of our detector decreases rapidly when the wavelength of 
light is shorter than 340 nm. Thus we could not resolve the subpeak at shorter wavelength of 
the double scattering spectra. In Fig. [3] (b) we compare the q value for the subpeak at longer 
wavelength to that of the secondary peak in the measured scattering spectrum. They exhibit 
similar dependence on 9. Therefore, the secondary peak in the measured light scattering 
spectra is mostly from the subpeak at longer wavelength of the double scattering spectra 
in Fig. m^b). As shown in Fig. [2](b), the measured secondary peak in the cross-polarized 
scattering spectrum is broadened and red-shifted compared to that in the co-polarized one. 
This is an evidence of triple scattering of light that causes stronger depolarization than 
double scattering. A qualitative analysis of three successive single scattering events reveals 
that the triple scattering spectrum is broader than that of double scattering and its center 
is in between those of single scattering and double scattering. 

Note that the secondary peak cannot be generated by high-order diffraction. For instance, 
the second-order diffraction gives ko = kj + qi + q2. It does not have such a constraint 
|ki + qi| = k, which applies to the double scattering due to k conservation in each of 
two independent single scattering events. Consequently the second-order diffraction peak 
is much broader in frequency. For example, in the backward direction ko = — kj, with the 
approximation I{q) ^ 6{q — qo), the value of k for the second-order diffraction ranges from 
and go- Hence, the second-order diffraction spectrum spreads from A = 280 nm to infinity, 
completely different from the measured secondary peak in Fig. [2](b). 

In summary, we demonstrate that double scattering in a short-range ordered nanostruc- 
tures produces additional peaks in the light scattering spectrum that contribute significantly 
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to structural colors. This phenomenon has no analog in either periodic structures or random 
structures. In a perfectly ordered structure, multiple diffraction or high-order diffraction 
does not create additional peaks, because a set of primary Bragg peaks would diffract into 
the same set of primary peaks. Vukusic et al. discovered new colors being generated by 
double reflection from the concaved multilayers of butterfly wings 16j. The normally inci- 
dent light, reflected from on one 45° inclined surface is directed across the concavity to the 
opposite orthogonal surface from where it returns parallel to the incident direction. This 
process strongly depends on the macroscopic geometry, and the direction of incident light. 
In contrast, the double scattering phenomenon we report here originates from the intrinsic 
nanostructure with short-range order. It is insensitive to the surface profile, and exists for 
any direction of incident light. Therefore, the creation of additional spectral feature by dou- 
ble scattering is a unique property of short-range ordered structures, that has been widely 
utilized by nature for structural coloration. 
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